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Abstract: Single Photon Detectors are integral to quantum optics and 
quantum information. Superconducting Nanowire based detectors exhibit 
new levels of performance, but have no accepted quantum optical model that 
is valid for multiple input photons. By performing Detector Tomography, 
we improve the recently proposed model [M.K. Akhlaghi and A.H. Majedi, 
IEEE Trans. Appl. Supercond. 19, 361 (2009)] and also investigate the 
manner in which these detectors respond nonlinearly to light, a valuable 
feature for some applications. We develop a device independent model for 
Single Photon Detectors that incorporates this nonlinearity. 
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1. Introduction 

Single Photon Detectors (SPD) have been central to our exploration of the fundamental limits of 
imaging fTl, interferometry ||2l, and communication f3|. They are used in astronomy lU, laser 
ranging [fSl, and biological imaging |6j. The SPD is a building block from which most more 
complicated optical detectors can be built. As examples, a photon number resolving detector 
can be built by multiplexing many SPDs Q, a phase sensitive detector by mixing the input light 
with an auxiliary laser beam |8|, and an entanglement (i.e. Bell-state) detector by coincidence 
counting between SPDs behind a polarizing beamsplitter |9|- Consequently, SPDs have been 
ubiquitous in quantum optics and quantum information experiments performed up to this day. 

Throughout, they have been described by a simple model. SPDs are binary detectors ( 'Click' , 
or 'No Click'), and thus any non-zero number of detected photons will result in the same re- 
sponse: a Click. This behavior is contained within the standard model for the SPD's positive- 
operator-valued measure (POVM): 

1(1 -A)"! «)(«!, (1) 

where tTq™ is the No Click operator. Pi is the quantum efficiency, | n) is an n-photon state, 
and the Click operator is tt^™ = 1 — tTq™. This simple model neglects the optical or electrical 
nonlinearities that likely exist, to some degree, in all SPDs. Nonlinear quantum optical detectors 
can prepare novel quantum optical states jTol, characterize the photon statistics of input states 
Q, and measure the temporal profile of optical wavepackets through intensity autocorrelation 
ifTTl . Advanced SPDs, with their sophisticated structures 1(121 . are good candidates to violate 
Eq. ([T]i through non-independent photoemissions and/or nonlinear gain mechanism. 

A good example is Superconducting Nanowire SPDs (SNSPD), which, compared to con- 
ventional SPDs, offer higher maximum count rate, lower timing jitter, broader spectral range, 
and superior detection efficiencies at telecom wavelengths lfT2l . These detectors consist of a 



meandering superconducting nanowire carrying a constant bias current close to their critical 
current. A single photon deposits enough heat to create a hotspot that is not superconducting. 
With the help of the bias current this can initiate a resistive bridge across the nanowire with 
some probability. The resulting voltage spike signals the detection of a single photon 11311141 . 
Unfortunately, apart from such descriptions and calculations, there is no associated full quan- 
tum optical model (i.e. valid for multiple incident photons) from which we could formulate a 
general POVM for a single element SNSPD. 

In an attempt to quantitatively describe the different observed experimental results within a 
single formalism, Akhlaghi and Majedi fl5\ proposed a semi-empirical model (AM model). In 
this model, the absorption of two or more photons at one point in the wire could result in coop- 
erative resistive bridge formation. Hence, one might expect a nonlinear response in SNSPDs. 
Moreover, it might dominate the linear response at low bias currents. This is in agreement with 
the early experimental results on characterization of SNSPDs 1 16 1. The model has been shown 
to work well with input light in a coherent state. 

To further test this model we perform tomography on a SNSPD. Detector Tomography ifTTl 
[TSl is an agnostic procedure to determine the POVM of a detector. The detector is treated as 
a Black Box in that we do not need to know its mechanism or make ancillary assumptions 
about it. Just recently demonstrated in principle. Detector Tomography has yet to be applied to 
a detector without an accepted model for its POVM. 

2. Detector Tomography 

To perform tomography we create a coherent state with a 1310nm laser that produces 200 ps 
pulses at 100 kHz. We use a programmable attenuator to set the coherent state amplitude, a. 
We sequentially send Rt = 10^ copies of this state into the SNSPD system and record the 
number of detector Clicks, Ri. This is repeated for a set of D states {| a,)}, increasing a from 
«() = until the detector response is unchanging at «£>. i.e. dRi /da = 0. The estimated CUck 
probability, R\/Rt, is equal to the Q-function of the Click POVM operator, Q{a). This, in 
itself, completely characterizes the detector response. 

Our SNSPD consists of 4 nm thick by 120 nm wide Niobium Nitride meandering wire on 
sapphire. The active area is square with a fill factor of 60%. The SNSPD is coupled to an optical 
fiber and cooled down to 4.2K by exploiting a cryogenic setup (described in |fT9l ). Fig. [Ha) 
shows a schematic of the measurement setup. The bias current is applied using a bias-Tee. The 
weak response of the SNSPD after amplification is split between an oscilloscope and a counter 
(see II20I for more detail). We use a 20ns counting gate, triggered by the laser, to reduce the 
dark count contribution from the times between input pulses. 

We perform Detector Tomography (i.e. measure Q-function) at three bias currents, 25, 20, 
and 16/iA. In Fig. [Ub), we plot the measured response for each of these (blue circles). We 
expect a standard linear SPD response at 25/iA since this is the normal operation mode. Using 

Eq. ([T]i, one expects a Q-function of the form Q{a) = 1 — exp ^— |ap^ . We estimate Pi 

using a single data point at Q{a) = R\/Rt =0.1. Indeed, using this Pi, the resulting predicted 
response (red line) agrees well with the measured response. Repeating this analysis for 20 
and 16/iA, we find that the estimated Pis, and thus, quantum efficiencies decrease as the bias 
current decreases. More significantly, the disagreement between the shape of the predicted and 
measured Click Probability distributions is substantial. Evidently, a SNSPD quickly becomes 
nonlinear as the bias current is lowered. 

A general limitation of Q-functions is that they cannot be used to calculate a detector's re- 
sponse to arbitrary states of light ||2 1 1 . Instead, one typically uses the POVM operators repre- 
sented in the photon number basis. We find the SNSPD POVM by fitting it to the D measured 
statistics R\/Rt according to the methods given in f22\. Since the detector lacks phase sensi- 
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Fig. 1. (a) Schematics of the experimental setup, (b) Click Probability at different bias cur- 
rents. The critical current was measured to be 26.0±0.5/xA. The blue circles are measured 
Click probabilities (not all points are shown). The red lines are calculated using the linear 
model in Eq. ([T) by calculating the efficiency parameter Pi when the click probability is 
equal to 0.1. The black dotted lines are from the nonlinear SPD Model in Eq. l|5j. {P,,} 
equals {7.30e-4, 2.49e-3}, {9.72e-6, 7.15e-5, 8.14e-9}, {0, 7.33e-8, 2.87e-10, 2.81e-14} 
for 25, 20, and l6flA, respectively. 



tivity, the off-diagonal elements of the Click operator are zero. We thus represent it as a vector 
n (with dimensions x 1). We truncate n at a number state — 1 that is sufficiently high that 
U{N) w 1. 

One difference with ||22j is that these particular SNSPDs have low system quantum efficiency 
(0.2% or less at 1310nm) and thus the maximum photon number, N ~ I, required to span their 
response was large. Instead of using large matrices in the fitting, we scaled the inputs {| a,)} 
by a factor <C 1. For each bias current, k is chosen so that the Click Probability is 95% at an 
average photon number (n) = 30. This scaled data is shown in Fig.|2a) (black circles). We plot 
n determined from it in Fig.|2jb) (blue circles). Using this n, in Fig.|2ja) (blue line) we plot the 
predicted detector response to coherent input states. This fits the scaled data well, confirming 
the fitting procedure. 

If scaled by too large a factor the POVM will, at the very least, be unrepresentative of the 
detector and, at worst, be physically impossible. Thus we test the estimated POVM for validity 
by adding k loss back into the POVM (i.e. Tl„„scaied = L^^^saiied, where L is the binomial dis- 
tribution matrix described in ll22l ) and predicting the detector response to the original unsealed 
inputs, {| a,)}. For all the bias currents, the difference between the predicted click probability 
and the raw data (i.e. \R\/Rt — Tr[n,„„f.a/g^/ |a,) (o!,|]|) is less than 0.15% on average and has a 
maximum of 1.4%. This indicates that we have accurately estimated the SNSPD POVM using 
the scaling technique. 
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Fig. 2. (a) The Q-function (i.e. click probability for coherent state inputs) of the scaled 
detector at different bias currents. The raw scaled data (black circles) agrees well with the 
tomographic POVM (blue line) and the AM Model (red dashed line), (b) The correspond- 
ing Click POVM operator. The operator found from tomography (blue circles) disagrees 
with that from the AM Model (red dashed line) but agrees with that from the nonlinear 
SPD Model (black dotted line). In the latter, {P„} equals {7.29e-4, 9.95e-2}, {1.08e-5, 
4.76e-2, 3.74e-3, 1.13e-4}, {0, 1.97e-4, 2.01e-3, 4.87e-4, 5.07e-5} for 25, 20, and 16/iA, 
respectively. 

3. Generalized Model 

The lack of ancillary assumptions and models in Detector Tomography make it general and 
objective. At the same time, it provides a surplus of information (i.e. O (N) parameters) that 
can be difficult to interpret. Tomography can hardly replace the natural ease and intuition that 
is associated with a model. A good model candidate is the AM model. Indeed, plotting the AM 
model in Fig.|2a) (red dashed line) shows that it can be fit to the scaled data. 

To find an analytic expression for the SNSPD POVM in the number state basis we repeated 
the derivation of the AM model given in fT5\ with one change: We derived the detector response 
to a Fock state rather than a coherent state input. The parameter values in the model should 
remain the same in either case. Thus, using the parameter values from the AM model curves 
in Fig. 12 a) we plot the AM POVM in Fig. IHb) (red dashed line). Substantial disagreement 
between the tomographic POVM and the AM POVM suggests there is an error. 

Indeed, further investigation revealed that in Eqs. (4) and (5) in ifTSl the average photon 
number is used in place of the actual photon distribution. The resulting error is only revealed 
in the comparison of analytical and tomographic POVMs. We fixed this error and found that 
the resulting analytic POVM has a form that is applicable to general SPD detectors. Leaving 
behind the SNSPD specific derivation in ifTSll . we now introduce a broad detector model that is 
able to describe nonlinear SPDs. We do this by generalizing Eq. ([T]) to multiphoton detection. 

Consider a binary detector that is only sensitive to n number of photons; any less and the 
detector responds with No Click, any more and it still only outputs one Click. This is the n- 
photon generalization of the SPD and, hence, we call it an n-photon detector (NPD). If m > n 



photons impinge on the detector, there are m choose n ways for those m photons to trigger the 
NPD. Consequently the generalization of Eq. ([T]i is: 



<™= E (2) 

m=0 

where tTq'™ is the No Click operator, P„ is the n-photon detection efficiency and ('") is the 
binomial coefficient (= for n > m, = 1 for n = 0). This generalization works even for a zero 
photon detector. We can identify Pq as what is commonly called the 'dark count probability'. 

A nonlinear SPD can be modeled as concurrent NPDs. As shown in Fig. |3ja), the model 
consists of a logical OR between M NPDs, where M represents the maximum number of mech- 
anisms that should be present in order to describe the response of the detector before it saturates 
at high intensities. The associated POVM operators are: 

oo M-l 

<^=i - L n(i-^")^"^NH, (3) 

m=0 /j=0 

for CUck and Kq^ = 1 - Trf^ for No Click. The CUck probability of the nonlinear SPD in 
response to a coherent state | a) is: 

Pcuckioc) = 1 - ^ e-l«l 1^ n (1 -P")^"^- (4) 

m=0 n=0 

This can be rewritten in matrix form as: 



C = FYl = F{E-exp{GH)), (5) 

where C (with dimensions Dx 1), includes the D measured statistics, Ri/Rj-; Yl (N x I) 
includes the diagonal elements of K^^; F contains the D coherent state probes, Fij ~ 

|a,p-'exp ^— |a,p^ £ is a matrix of ones; G is a matrix of binomial coefficients such that 

Gij = (jZi); ^nd H (M X 1) is an unknown matrix which includes the unknown set {Pn}, 
=ln(l-P„=,_i). 
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Fig. 3. (a) The Nonlinear Single Photon Detector model. Each element represents an n- 
Photon Detector. A broadly applicable model is created by logically ORing these elements, 
(b) {P„} at 16^A found under different scaling factors, k. 0, Q, □, y, A, and -{f represent 
P„=0 to Pn=5 respectively. The dotted lines have slopes equal to n. The plot shows a linear 
optical loss of rj scales P„ — >• ri"P„. 



We estimate H by solving the following constrained nonlinear multivariable optimization 
problem: 

C-FE+Fqx^{GH) 



(6) 



C 

subject to // < 0. The second norm of a matrix is defined as ||A||2 = I Li,/ V^i.n ) ■ Each 
element of the expression is normalized to C to give equal weighting to all the points. The 
constraint of the problem ensures the optimization leads to a physical result for the set 
We note the function exp(ax) is convex on R for any a G R fSSl, which also makes Eq. ^ 
convex. 

We first solve Eq. ^ for scaled input states. From the estimated {f„} we only keep those 
elements that change the minimum of Eq. ^ by more than 1% (25/iA: {fo, ^"1}; 20/iA: 
{Po, • • • , A}; 16;UA: {Pi, ... , P4}). These parameters classify the operation of the SNSPD, 
from a standard SPD at 25/iA to a composite of one, two, three, and four photon detectors at 
16M- 

From these {Pn} we calculate the nonlinear SPD Click POVM operator for the three bias 
currents in Fig.|2jb) (black dotted line). They agree with the tomography POVMs to within 1% 
for most of the elements of H. For 25jUA and 20/iA the maximum difference is 3% at n (n = 1 ) 
and 6% at n [n ~ 3) for 16/iA. This excludes the large error at n (« = 0), which we attribute to 
insufficient measured statistics at extremely small mean photon numbers. The Quantum Fidelity 
(see 1221 for a definition) of the model and tomography operators are above 99.8% for all 
three bias currents. Thus, the model successfully gives the POVM of the SNSPDs including 
their nonlinearity, but with a dramatic reduction in the number of parameters compared to 
tomography. 

We can attribute the SNSPD nonUnearity to the fact that at lower bias currents the absorption 
of multiple photons in close proximity and within a short period of time can inject enough 
energy to switch the wire to normal state more efficiently than independent absorptions. 

In Fig.[TJb), we plot the model's predicted response to the coherent state inputs, with param- 
eters from the unsealed data (black dotted line). The difference with the measured statistics is 
less than 0.34% on average (i.e the minimum found in Eq. (|6]l divided by the number of points). 

Unlike the standard SPD, where Pi — > 77 Pi under a preceding optical loss of 77, there is no 
analytic formulae for how {P„ } transform under loss. By inspection of the scaled and non-scaled 
model fits (see Fig.[3lb)), however, each element of {P„} that is significant approximately scales 
as Tj". So removing any linear optical inefficiency from a nonlinear photon counter makes it 
more nonlinear. We also note Pq is not dependent on the optical input as expected. 

4. Conclusion 

The almost identical outcome of Detector Tomography and the nonlinear SPD model confirms 
both of them are reliable. This is a good example of how Detector Tomography is particu- 
larly useful for characterizing detectors outside their normal operating regime, where there is 
no model for their operation. We expect the nonlinear SNSPD model to be useful for other 
nonlinear binary detectors such as two-photon absorbing Avalanche Photodiodes and Electron 
Multiplying CCDs (thresh-holded). It will also be useful for characterizing conventional SPDs 
and putting limits on their nonlinearity. 
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